function f = fun9(x)
f = [2*x(1)^2 + x(2)^2 - 48*x(1) - 40*x(2) + 304
    -x(1)^2 - 3*x(2)^2
    x(1) + 3*x(2) - 18
    -x(1)-x(2)
    x(1) + x(2) - 8];
